Effects of a dipolar field in the spin dynamics of a Fermi liquid 
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We study spin dynamics of a normal Fermi liquid taking into account the demagnetizing field 
produced by the spin system itself. Linear solutions of the spin dynamics equations in the form 
of standing spin waves in a finite volume of liquid are found. At almost all known experimental 
conditions the influence of demagnetizing field can be satisfactorily described by the first order of 
perturbation theory. We carried out perturbational calculations for two geometries of experimental 
cell — spherical and finite-cylindrical. We performed also exact numerical simulations of the spin 
wave spectra in a spherical cell at an arbitrary strength of the demagnetizing field. 

The obtained results are applied in particular to conditions of recent experiment (G.Vermeulen 
and A.Roni, Phys. Rev. Lett. 86, 248 (2001)) related to the problem of zero temperature transverse 
relaxation in a polarized Fermi liquid. We found that not taking into account demagnetizing field 
leads to negligible errors in the measured relaxation time, thus supporting the conclusion of the 
absence of zero temperature spin wave damping. 
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I. INTRODUCTION 

Spin dynamics of a strongly spin-polarized normal 
Fermi liquid still captures appreciable theoretical and ex- 
perimental interest. Among the main questions here is 
whether the transverse (i.e. in a direction perpendicu- 
lar to the external field) magnetization excitations are 
damped at zero temperature. 

Polarizing a Fermi liquid creates a gap ~ h^/M/xn 
between the two Fermi energies for spins up and down. 
Here 7 is the gyromagnetic ratio of 3 He nuclei, Xn is 
the susceptibility of the liquid. Meyerovich has pointed 
out that the existence of the gap leads to a non- 
conventional temperature dependence of the transverse 
relaxation time, oc (T 2 + where T a is of the or- 

der of the gap, and therefore to a damping of the trans- 
verse excitations even at zero temperature. This idea has 
been pursued in several theoretical papers 0, JiJ, [Q and 
was contested recently by Fomin ^ , who argued that the 
conclusion of existence of zero-temperature attenuation is 
drawn from the wrong premises about the ground state of 
a polarized Fermi- liquid, viz., from treating the quasipar- 
ticles between two Fermi levels as excitations. Whereas 
as long as the polarization of the liquid is conserved, these 
particles should be considered as an inalienable part of 
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the ground state. 

Curiously, a similar discussion arose in mesoscopic 
physics where a poorly argued concept of finite dephas- 
ing time at T = has been developed and contested 

S El- 

The results of the measurements of the spin diffusion 
coefficients by spin echo experiments in pure 3 He ||, 
and in solutions of 3 He in 4 He J8), J9) with a spin po- 
larization of a few percent revealed a finite value of T a 
in a qualitative agreement with the zero-attenuation con- 
cept. But the observed T a were several times greater than 
theoretical estimations in ||, ||, On the contrary, 
the recent measurements of linear spin wave damping in 
dilute 3 He at even higher polarizations jlO| are in agree- 
ment with Fomin's theory (T a = 0), although the upper 
limit for a finite T a set by this experiment (due to the 
error bars) does not allow to rule out completely the ex- 
isting theory of zero temperature spin wave damping. 

A coherent theory of strongly polarized Fermi liquids 
based on a properly defined ground state is lacking. On 
the other hand a proper interpretation of the experimen- 
tal data for a strongly polarized liquid is itself non-trivial. 

The point is that the magnetic field acting on the spins 
of a liquid is conventionally supposed to equal the ex- 
ternal field H e . In reality the field inside a specimen is 
well-known to differ from H e due to the shape-dependent 
demagnetizing field proportional to the magnetization. 
An oscillating magnetization thus acts back on itself via 
the demagnetizing (or dipolar) field. This phenomenon 
manifests itself as magnetostatic waves in ferrimagnets 
theoretically described by Kittel ||l4| and Walker |fl5|| . 
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Walker originally showed that if the magnetization M 
is supposed to obey the Larmor precession around the 
internal magnetic field LP 

(d t + 7 H*x)M = 0, (1) 

then because M and H? are related through Maxwell 
equations they can self-consistently oscillate only with 
certain frequencies localized in the range ~ 2ir~/M near 
the Larmor frequency. 

In an interacting Fermi liquid (|l|) should be replaced 
with the Leggett system of equations. In the linear ap- 
proximation, solutions of this system are standing spin 
wave modes with the widths proportional to the trans- 
verse relaxation time. So study of the behavior of the 
widths of the modes with temperature is one of the pos- 
sible ways to detect the zero-temperature transverse at- 
tenuation. 

For weakly-polarized liquids the effects of the demagne- 
tizing field can be discarded. But it is preferable to have 
strong polarizations in order to increase the predicted 
temperatures of the attenuation onset. So the effect of 
polarization has to be taken into consideration for the 
proper interpretation of the spectra. On the other hand 
at sufficiently strong polarizations one can expect con- 
siderable changes in the Leggett description of the Fermi 
liquid spin dynamics. In particular, due to the presence 
of two Fermi surfaces a double set of the Fermi-liquid 
parameters must enter the theory. We will use Leggett 
equations assuming that they are still valid when the 
share of polarized nuclei of the liquid does not surpass 
~10 %. 

To include dipolar field we have chosen to write out 
the dipolar part of the internal field B? explicitly as 
an integral of the magnetization, this integral being a 
general solution of Maxwell differential equations with 
Maxwell boundary conditions. Thus we work with a 
closed integro-differential equation directly on the mag- 
netization. 

A short review of other possible methods is done in the 
Discussions section. 

We start directly from the generalized Leggett equa- 
tions and study the effects of the demagnetizing field co- 
herently. We specialize to the case of linear spin waves, 
setting spin echo experiments aside. 

A full-blown study including numerical calculations is 
done only for a spherical shape. For a finite cylinder 
as well as for a sphere we also calculated the changes to 
the spectra by the demagnetizing field using perturbation 
theory. 

Our results show that the demagnetizing field intro- 
duces small corrections (about 4%) to the value of low 
temperature transversal relaxation rate experimentally 
determined by G.Vermeulen and A.Roni [[To). Hence, the 
main conclusion of Ref. about the absence of zero- 
temperature spin wave damping is supported. 

The paper is organized as follows. Sections II and III 
form the basis needed to comprehend the authors' point 
of view. In Sec. II we show how to include the dipolar 



field in the standard Leggett equations and to how to 
linearize the result to obtain an equation for spin waves 
subject to both exchange and demagnetizing fields. 

In Sec. Ill as the simplest application of the theory 
developed we find corrections to the spin wave spectra in 
a finite cylinder in the first order of perturbation theory. 
At the end of the Sec. Ill we find dipolar limitations 
on the correct determination of the transverse relaxation 
time from the conventional interpretation of the spectra. 

The next Section IV contains similar first order per- 
turbational estimation of the dipolar-field corrections to 
the spin waves spectra in a spherical container. 

The results for a sphere are compared to numerical 
simulations carried out in Section V, where we also cal- 
culate spin wave spectra in the regimes of intermediate 
and strong demagnetizing fields. 

In the last Section VI we discuss conclusions. 

II. STATEMENT OF THE PROBLEM 

A. Basic equations 

Spin dynamics of a Fermi liquid is described by the 
Leggett coupled system of two partial differential equa- 
tions on the local magnetization M(r, t) and its current 
J<(r,t) © 

(dt+jBx) M +diJi = 0, (2) 

2 

(dt+jBx) J, + ^(M-Mo) (3) 
+ k— MxJ, = -- • 

Xn Tl 

Here B is the flux density inside the sample j2lj], w 2 the 
renormalized Fermi velocity w 2 = V F (1 + F§ )(1 + Ff/3) 
and ti the renormalized relaxation time t\ = r/(l + 
Ff/3). Fq and Ff are the coefficients of expansion of 
the antisymmetric part of the Fermi-liquid interaction in 
the spherical harmonics. The equilibrium magnetization 
is 

Mo = Xn H*, (4) 

where FT = B — 47rM is the internal field. 

Spin dynamics equations of a Fermi liquid reduce to the 
form (H) , (|^) in either of the regimes — collisionless C ^> 
1 or hydrodynamic C <C 1, where the regime parameter 

C = K(7iF)(M/Mo)Ti. (5) 

The factor M/Mq accounts for the possibility for the po- 
larization M be hig her than the equilibrium value Mq — 
in experiments Jl(|, |l8| M/Mq varied from 1 to 5. 

The condition of applicability of the Leggett equations 
|pp| is that the characteristic scale of spatial inhomogene- 
ity £ be greater than the quasiparticle mean free path vft 
or the magnetic length vpr/C which one is the shorter, 

£ > min{v F T, v F t/C} . (6) 
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In the collisionless regime (C > 1) the magnetic length 
vpr/C is the shorter and thus should be smaller than £. 
While in the hydrodynamic regime (Cc 1) the spatial 
scale £ should exceed the mean free path vft. 

Eq. (||) contains the torque due to the local molecular 
field «M(r, t)/xn acting on current Jj(r, t). The combi- 
nation of the Fermi-liquid constants Fq and F( L 
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(7) 



It 



measures the strength of the exchange interaction 
vanishes when turning the exchange off. 

Leggett originally |2(J considered the case of a weakly- 
polarized sample, wherein B = H e — the external mag- 
netic field at infinity. Generally, the relation between B 
and H e is to be determined from the conventional bound- 
ary value problem of solving the magnetostatic equations 
in a non-conducting medium 



d x H 



0. 



9B = 



(8) 



with Maxwell boundary conditions of the continuity of 
B n and H t at the boundary of the sample and of H — > H e 
at infinity. In (pt) 



H = B 4ttM. 

A general formal solution of 
propriate boundary conditions is |22| 




H H c + Hdip, 



(9) 

satisfying the ap- 
(10) 



where 



H dlpW = 8 ( 8 /M^^ 



(11) 



is called the dipolar field. It is straightforward to verify 
that Eqs. ( fic| ) and @ are indeed the solution in the whole 
space with the help of 



d 2 \v 



(12) 



M in its turn has to be found from the Leggett equa- 
tions (||), (0). Therefore, the closed system of integro- 
differential Eqs. (g), (§), (|)-(|l|) completely describes 
normal Fermi liquid electrodynamics with the effects of 
both inhomogeneity and the demagnetizing field taken 
into account. 



Inside the sample the field H from (10) is called B? 
— the internal magnetic field. The difference between 
external field at infinity H e and H l is usually denoted as 

H* — H e = — 47rn[M], (13) 

where the (tensor) operator n, acting by the rule 

(14) 



m] = -la(of v 



is called demagnetizing operator. Let us agree to denote 
an operator with a hat over a letter and a tensor with 
an underlined letter. The demagnetizing tensor opera- 
tor n is specific to the shape of the sample, over volume 
of which the integration in (|l4]) is taken, and generally 
is coordinate-dependent. It only reduces to the tensor 
of constant demagnetizing coefficients n when acting on 
spatially homogeneous distributions in ellipsoidal sam- 
ples (including limiting cases of a slab and an infinite 
cylinder). 



B. Static magnetization distribution 

In a linear spin wave the magnetization rotates with a 
small amplitude m around its (large) stable value M in 
a static external magnetic field H e = H e z. 

Consider a general static (<9 t M = 0, J; = 0) case. In 
order for J,- = there should be 



Si(M-M o )=0 



(15) 



and then for <9 t M = the magnetization should be locally 
directed along H' 



(16) 



M(r)= Xn A(r)H 4 (r). 



The function A(r) must be such that M(r) satisfies (|l5|). 
This imposes restrictions on the spatial dependence of 
A(r), but leaves A(0) arbitrary. 

So there exists continuum of static non-equilibrium 
magnetization distributions numbered by A = A(0). It 
is this A which is represented as M/Mq in Table I. 

In equilibrium (Q) A(r) = 1. 

To find a form of a non-equilibrium static magnetiza- 
tion distribution, we shall use the smallness of Xn- For 
pure 3 He, the magnetic susceptibility is 



Xn = n VJVb/2(i + iff) ~ lo- 



(17) 



where A*o = m*kp /2n 2 fi 2 is the density of states on the 
Fermi surface. For 3 He- 4 He mixtures Xn is less, propor- 
tional to kF oc \fx, where x is the concentration of 3 He 
atoms in the mixture. 

Substituting ((TJ) into @ yields 

H*(r) = H e (r) - 4^ X „n [A(r)H e (r)] + 0( X 2 n ). (18) 

For the reasons that will become clear below the exter- 
nal magnetic field is taken almost constant, with a small 
gradient along its direction 



H e (r) =li e {l + zVu l /lo l ). 



(19) 



The presence of the field gradients in the perpendicular 
directions, necessary for the fulfillment of the condition 
dH e = 0, is inessential for the following discourse. 

From @ it then follows that VA = (1 - A)Wuj l /lu l , 
i.e. the spatial inhomogeneity of A(r) = A + zVA has 
the same smallness. 
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Leaving in ( |l8| ) only the first order terms in either Xr, 

or zVlol/^l, we have 



Here 



7H l = {ujl + z\7ujl)z — ujmu[z\. 



lu l = 7-fF 



is the Larmor frequency and 



(20) 



(21) 



(22) 



is the frequency corresponding to magnetization. n[z] 
is generally a coordinate dependent vector. For ellip- 
soidal samples it reduces to nz. If one of the principal 
axes of the ellipsoid (which are also the principal axes of 
the demagnetizing coefficients tensor n) coincides with 
z, we have nz = zn^ z \ where nS z > is the z-th demag- 



netizing coefficient. E.g., for a sphere n 



0) 



for a 



plane-parallel slab with the edges perpendicular to z the 
coefficient n^ = 1, and for an infinite circular cylinder 
with the generatrix parallel to z the coefficient rS^ = 0. 

The small ratio 2RWujl/^l, where 2R is the sample 
size, together with 47rx n ^i(0) ~ 10 -6 are the two small 
parameters in the problem. Conventionally the dipolar 
field is not taken into account and the second parameter 
is considered negligibly small. This is no longer justified 
for recent experimental conditions as is seen from Table 
I, where the ratio of the second parameter to the first 

w m /2RVuj l (23) 

is represented for various experimental conditions. 

C. Linearized equations of motion 

To obtain the linearized form of the equations of mo- 
tion (||), (||) we expand all the macroscopic quantities 
near their stationary values: 



H e = H e z- 
M = M 



m. 



H 1 
J, 



H* 



(24) 



where the static value of H l is that from ( p0| ) . The radio 
frequency field h e plays the role of a driving force for the 
spin system response m. We denote h 1 = h e — 47m[m]. 

The static M in (|24|) is the result of the substitution 
of (|2C|) into (Hi 



jM/xn = uj l {A + zWoj l )z - uj M AH e n[z}. (25) 
However, we can retain only the greatest term in M 



M = Xn An e 



(26) 



when linearizing Eqs. (g), (||). In (||) this is simply due 
to the fact that the terms of M to be omitted are of the 
order 0(x n zVuj L /uj L ) and 0(x„)- In © M is divided 
by Xn and the justification is lengthier. We will suppose 



and discuss why only the main term in M should be 
left after the derivation below. 

In practice one usually is interested in movements qua- 
sistationary in the Larmor frame. To a first approxima- 
tion one supposes that jj(r, t) is stationary, i.e. precesses 
with the frequency 7B: (<9 t +7Bx )jj = 0. Then resolving 
Eq. (0) with respect to ji with M from (Eq) gives 



Ji = 



D 



1 



^2 [dun -Czx 9,m + C 2 z{zdijn)] , (27) 

\v%t{1 + 



where the diffusion coefficient Dq 
Fg) and 



1 2 
3" 



C = KTi-yM/xn = ktxAlol 

is another expression for the regime parameter 
One then plugs the divergence of (|27| ) into 1 
divergence of ji has the order 



diji 



D d 2 m 



(A)/£ 2 )m, 



(28) 



The 



(29) 



where £ ~ ^/ Dq / Vwi is the characteristic scale ( |49|) . So 
ddi ~ (^Vwi)m has already the smallness zVljl/ujl. If 
we had accounted for the higher order terms in M than 
( p6| ) when calculating the current, these terms would 
have entered ji through the regime parameter C, and 
after multiplication with Dod 2 m would have produced 
terms of the order 0(x n zV 'u> 1 / 'to 'L)m and 0(Xn )m. That 
is why we were allowed to substitute simply ( p6| ) in (|^). 

In the linear approximation m(r) is in each point per- 
pendicular to the static M(r) if the absolute value of the 
magnetization is conserved. Exp. ( p5| ) shows that apart 
from the major x— and y— components m also has a 
minor m z ~ Xrjn x - This component also may be seen to 
give higher order terms and is therefore negligible. 

We will thus consider m_Lz. The last term on the 
right-hand side of ( |2^ ) then vanishes, and substituting 
(|3) into (§) one gets 

(d t - §<9 2 ) m + z x (-/H l z m + Dd 2 m - 7 A/h*) = 0. 

(30) 

Here we left only the z— component of B? because Hj_, 
which multiplies vectorially only by m z z, produces terms 

°(xl). 

In ( 130] ) we introduced the effective spin diffusion coef- 
ficient 

D = D C/(1 + C 2 ). (31) 

In the strong (C ^> 1) collisionless regime D w 
w 2 /SkujlA is temperature independent. 

Suppose m, h depend on time as monochromatic waves 



m, h 



oc e 



-iwt 



Then written out in components the equa- 



tion of motion for m becomes 

.~ (m x \ f-m v 
-iuj\ + UJ L \ 



um f —h 

47T 



(32) 



where we designated by Si, and Q respectively the oper- 
ators 

Q L =jHl + Dd 2 , Q = uj-i^d 2 . (33) 
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Multiplying by 



1 

-1 



results in the second order inhomogeneous partial differ- 
ential equation £ _1 m = h^_, or 



X 1 m + 47rn[m] = h^_. 



Here x 1 is the inverse susceptibility tensor 



V 



47T 



iu) u)l 



(34) 



(35) 



To close the boundary value problem one must impose 
some appropriate boundary conditions on m. Supposed 
that the container is made from a non-magnetic material, 
there is no magnetization current into the walls, and we 
get boundary condition in the form 

fiidim\ d = 0. (36) 

Here fn is a unit normal to the wall. 

In the normal variables m± = m x ± im y Eq. (|34|) has 
equivalent form 

(Q L T w)m± 



mm 

47T 



is the zz-component of the demagnetizing tensor (|i4|). 
And n zz [l] is a scalar function on coordinates, which ac- 
tually coincides with the z-th demagnetizing coefficient 
n^ 2 ) for ellipsoids with one of the principal axes parallel 
to z. cj/,(r) denotes ojl + zVwj,. 

To conclude, we have derived equations of motion ( |34l ) 
or ( |38| ) for small deviations of magnetization from static 
values. The magnetization as a function of external rf 
field is a response of the system on a particular radio fre- 
quency to. The full form of the equations of motion ( |3^ ) 
is unnecessarily complicated because it contains superflu- 
ous information on the dynamics of the counter-rotating 
component m~ of the magnetization. An example of so- 
lution of the full Eq. (|34j) for an infinite medium in a 
uniform magnetic field is analyzed in Appendix |a|. 

The rest of the paper deals with ( |38| ) . 

The response of the liquid is detected through changes 
in the impedance of the NMR coil, which are proportional 
to (see Appendix [§]) 



x = (h e+ \g^\h e+ ), 

where is the Green operator 

|m+) = {oj M /^)Q L o\h e+ ) 



(41) 



(42) 



Green operator may be expanded into an infinite sum 
in the eigenfrequencies u) a of the homogeneous equation 
corresponding to (|38|) 



whence we estimate 

?7i_ ~ (tu — (jJi^m+^ujL -C m + . 

We see that in the vicinity of the Larmor frequency, 
when w ~ wl, the counter- rotating component to_ may 
be neglected with respect to the co-rotating m+. Then 
(U) simplifies to a single linear inhomogeneous integro- 
differential equation 



(H-u ) )m + (v) = ^hl(v) 



(38) 



with a generally non-Hermitian Hamiltonian operator 



Hm + = D(l+ ±)d 2 m++ 1 H l z {Y)m + 



8tt 



d + d- 



(39) 



We may use the equivalence d+d- = d 2 — d%, (|2^) and 
the property (|l2|) to rewrite the Hamiltonian as 



Hm+ = D (1 + i) d 2 m + + uj L (r)m+ 



+ 



((l-2n^[l])m + -n Z2 [m+]). (40) 



Here the integro-differential operator 

Z(r') 



H\a) = lu\o) 
with the boundary condition 



= 0. 



(43) 



(44) 



For a Hermitian Hamiltonian the eigenfrequencies ui a 
are real and the expansion is 



\a)(a\ 



to 



(45) 



The absorption spectrum in the Hermitian case consists 
of a series of <5-peaks at w = w a . Indeed, writing real uj a 
in ( |45| ) as u> a + iO we see that the imaginary (absorption) 
part of @ 

is a weighted sum of (5-functions 



■El 



(a\h e+ ) \ 6(u - uj a ). 



(46) 



In the general case of a non-Hermitian Hamiltonian the 
expansion ( ^5| ) should be revised. We postpone the ap- 
propriate discussion until Sec. V. Here it is enough to say 
that the numerators in the series ( p5| ) remain the same 
in the general case, but the eigenfrequencies Lu a become 
complex, meaning that in general spectrum consists of 
Lorentzians. 

The following important conclusion drawn on the basis 
of (45) also holds in the case of a non-Hermitian Hamilto- 
nian. In a homogeneous (Vu>l = 0) external static fields 



H e and for ellipsoidal samples the Hamiltonian (40) has 
uniform solutions, so-called Kittel modes iQ with the 
frequencies 



lu = uj l + — (I - 3n zz ) 



(47) 



where n zz is the zz-component of the demagnetizing co- 
efficient tensor n. 

For customary sample sizes the rf field h e may be con- 
sidered spatially uniform. Then from ([46]) it follows that 
it is impossible to excite a non-uniform mode by a homo- 
geneous rf field h e , because then (a\h e+ ) oc (a\0), where 
|0) is the Kittel mode, and different modes are mutually 
orthogonal (a|0) = S a Q. 

We conclude that in order to couple to non-uniform 
eigenmodes the external static magnetic field should be 
inhomogeneous (see fll9|)) so that there would not exist a 
uniform eigenmode. 



III. FINITE-CYLINDRICAL CELL 

Study of the eigenstates of the Hamiltonian ( |40| ) in 
general is possible only numerically. This has already 
been done in Ref. pi, [ |l7| neglecting the contribution 
of the dipolar field. The former work dealt with ( |4C| ) in 
rectangular boxes, while the latter — in spherical con- 
tainers. 

In this paper we study the dipolar field effects due to 
the third term in ( |4"(i| ) numerically in Section V. At the 
same time, in the case of ujm = the problem allows 
explicit analytical solution for some typical experimental 
conditions. Such solutions are of undeniable interest — 
with them in hand we may use perturbation theory to 
calculate corrections to modes in the first order in u>m- 
So the two following sections are dedicated to such solu- 
tions and to the calculations of perturbational corrections 
respectively in the geometries of a finite cylinder and a 
sphere. 

We start the consideration of finite cylinders from ide- 
alized one-dimensional geometry of a plane-parallel slab. 
Next we calculate the first order perturbations to the 
modes frequencies due to the finiteness of the cylinder. 



A. Slab 

In the absence of dissipation, when C _1 = 0, the ef- 
fective diffusion coefficient D(l + iC~~ r ) is real and thus 
the Hamiltonian Ti. ( |40| ) is Hermitian. 

In the slab geometry the solution should be sought in 
the form 

m+(r) cx e ?kri m+(z), 

where rj_ is the coordinate vector in the plane perpen- 
dicular to z. Nevertheless, as is clear from (f4l|), only the 
solutions with k = contribute to the observation signal 
if the rf field h e is homogeneous. 



The eigenfunctions m^{z) = (z\a) then are the com- 
binations of the two Airy functions 

(z\a) = AAi (arg) + BBi (arg) , (48) 

where arg = (u a — lol — zVcjl)/£Vojl, and 

(49) 



is the characteristic wavelength. Its sign depends on the 
relative sign of D and "S/ujl and thus on the sign of k. In 
3 He and in 3 He- 4 He solutions with a concentration x > 
3.5%, K is positive. We consider £ > for definiteness. 

The boundary conditions ( fl4|) on the two plane bound- 
aries d z |a)| 2=0 L = determine the eigenfrequencies ui a 
and the ratio of the coefficients BJA. The remaining coef- 
ficient A is determined from the normalization condition 
(a\a) = 1. 

When L ^> £ the influence of the lower wall of the 
container is negligible and the modes are localized near 
the upper wall and decay exponentially into the bulk on 
distances ~ £. Then the eigenfunctions are just the Airy 
functions of the first kind — Ai and (48) becomes 



,) = AAi ( — : -+■ 



where a' n < is the n-th zero of the derivative of the 
Airy function Ai': a[ w -1.02, a' 2 w -3.25, a' 3 w -4.82, 
etc. The eigenfrequencies are 



ut + LVw L + a' n £Vw L . (50) 



Inclusion of dissipation C _1 ^ makes the diffusion 
coefficient complex. The Hamiltonian ( fiol) then becomes 
non-Hermitian. The complete analysis of the spectra of a 
non-Hermitian Hamiltonian is possible only in the frame- 
work of the general formalism to be developed in Sec. V. 

However for the moment it is sufficient to make the fol- 
lowing statement. In the presence of dissipation (C _1 7^ 
0) £ in ( f49| ) becomes complex 



f-tfl + i/C-) 1 ^ 



(51) 



and so do arguments of the eigenfunctions ( [48| ) . Eigen- 
frequencies (|50| ) also acquire imaginary parts due to £ 
entering the expression. 

Thus the complex eigenfrequencies in the presence of 
dissipation can be easily obtained from the real ones in 
the absence by the substitution (||l]) . This statement ap- 
plies not only to ( |50| ) but to any spectrum of the Hamil- 
tonian ©. 

Let us now consider the effects of the dipolar field. 
We may utilize the results of Appendix |A] since a slab 
is infinite in the direction perpendicular to z and modes 
depend on only z. In such conditions the demagnetizing 
field is local 



''slab 



(52) 
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It is then obvious from ( j40| ) that the dipolar field pro- 
duces no effect on the spin wave spectrum ( |50| ) apart from 
a uniform shift by 



1 



ww(l-3n^ b ) = -wm 



(53) 



Such a shift does not distort the spectrum — it does 
not change neither the mutual positions of the modes nor 
their widths, from which the characteristics of the liquid 
are derived. 

We rather aim at finding those distortions, so we pro- 
ceed to a more relevant shape of a cylinder of a finite 
radius, which as well as all finite shapes as we will see 
does give such distortions. 



B. Finite cylinder 

Consider a finite cylinder with the base radius R » £ 
and height L £ and a generatrix parallel to z. We 
find the influence of the finiteness of a specimen on the 
magnitude of dipolar corrections to the spectrum in the 
first order of perturbation theory. 

The first order perturbational corrections to the modes 
frequencies are the averages of the perturbation (dipolar) 
operator in the given eigenstate ip a ( T ) 

*dip«a = ^(l - 2 J \Mr)\ 2 n zz [l]d 3 r 

C(r)aj^(r)]d 3 rY (54) 



In a finite cylinder the dipolar-free eigenfunction |a) 
satisfying boundary conditions (44), written in cylindri- 
cal coordinates z, p, ip is 



-"n z n m £ 



Ai(^+a' )J m (%t p ) e 



(55) 



where n Zl n p = 0, 1, 2, oo are respectively the longitu- 
dinal and radial quantum numbers and m = — oo, +oo 
is the azimuthal quantum number. Q' n m is the (n p +l)- 
th zero of the derivative J' m of the Bessel function J m , 
Cn z n p m are the normalization coefficients. 

From the general formula ( |il| ) it is not hard to see that 
only the modes with n p = m = couple to the homo- 
geneous rf field. Indeed, |n z 00) is uniform in the plane 
perpendicular to z. Therefore, integrals like (n z n p m\h +e ) 
are proportional to (n p m [00) = <5„ p o<5mO- 



Calculating (Q) with ( p5|) yields (see Appendix [c| 



5 dip u, c Z Vl nd ° r = % ( -1 + ^ log ^ 



L , 8R 
— log — 

ttR eL 

8R 



(56) 



The first two terms not dependent on the mode num- 
ber n z describe uniform shift of the spectrum, and the 



last two terms dependent on n z through the numerical 
constants <J>„ = , ^/ n , , E nz and 9„ s , which are of the order 
of unity (see Table]n] in Appendix [c| , give the sought-for 
spectrum distortion. 

We see that for finite £/ R the spectrum undergoes dis- 
tortion proportional to the parameter 



2uj m 



7t£Vwl R 



log- 



(57) 



where ujm = 4irjM characterizes the magnetization den- 
sity, Vu>l is the gradient of the Larmor frequency, R is 
the radius of the cylinder base, L its height and £ is the 
wavelength ( p9| ) of an Airy-type standing spin wave. The 
quantity (Vwl gives the average distance between modes 
in the units of frequency. $ and \& are numbers of the 
order of unity. 

The calculations of the dipolar field effects in the first 
order of perturbation theory allow us to estimate the 
maximum error due to the demagnetizing field in the 
determination of the transverse relaxation time r from 
the spin wave spectra (see Appendix |cj) . 
^This error turns out to be of the order of the parameter 
i.e. for the experiment w 4.2% for M/M - 4. 
Thus, interpreting spectra according to the usual the- 
ory not taking into account the demagnetizing field in- 
duces an error in the derived value of the transverse re- 
laxation time of the order of the parameter (|S7l) . 

The term proportional to <& n , comes from the demag- 
netizing field produced by the rotating part m of mag- 
netization. While the term proportional to is due 
to the spatial inhomogeneity of the demagnetizing field 
— 47rn[M] produced by the initial static (homogeneous!) 
distribution of M ( p6| ) in a finite cylinder. 

The values of these two terms are plotted in units of 
lom as functions of the mode number in Fig. [j] for = 
0.015 and L = 2R. Apart from being greater, the term 
brought about by the spatial inhomogeneity of the static 
distribution of the dipolar field depends stronger on the 
mode number, thus resulting in bigger mutual shifts of 
the modes. So the main source of the spectrum distortion 
in a finite cylinder turns out to be the inhomogeneity of 
the static dipolar field. 

In ellipsoids, in particular in a sphere, the demagne- 
tizing field produced by the initial static distribution of 
M is homogeneous. And so there is only the distortion 
to the spectrum from the rotating part m of magneti- 
zation as we will see in the next section. This makes 
ellipsoidal shapes advantageous if the dipolar field effects 
are unfavorable. 

To conclude, we found the corrections to the spin wave 
modes in a finite cylinder due to a weak demagnetizing 
field in perturbation theory. These corrections consist in 
shifting the spectrum as a whole, changing the relative 
distances between the modes and in narrowing down the 
modes. The two last are of interest for us since they 
deform the spectrum. 

There are two contributions to the spectrum defor- 
mation — one from the static inhomogeneous demag- 
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FIG. 1: Comparison of the two parts in the parameter of dis- 
tortion of the spin-wave spectrum in a finite-cylindrical cell 
brought about by the demagnetizing field (the two last terms 
of in units of ujm for £JR = 0.015 and L = 27?). The 
first part (third term in (p6|)) (Q) is due to the demagne- 
tizing field produced by the rotating part m of magnetiza- 
tion. The second part (the last term in (|56|)) (□) results from 
the inhomogeneity of the dipolar field produced by the ini- 
tial static distribution M = XnAH e z of magnetization in a 
finite-cylindrical sample. 



netizing field and the other from the rotating part of the 
magnetization. The first contribution exists only in non- 
ellipsoidal samples, in which a homogeneous static mag- 
netization produces inhomogeneous demagnetizing field. 



IV. SPHERICAL CELL 

Though for a spherical container exact analytical so- 
lution in the absence of dipolar field turns to be impos- 
sible, one can obtain an explicit expression for several 
first modes in adiabatic approximation if the radius of 
the sphere R ^> £. 

The prerequisites of adiabatic approximation might 
be best understood if one exploits the analogy with the 
Schrodinger equation for a particle moving in an exter- 
nal field. If the movement in one direction is somehow 
more restricted than in the others (geometrically or by 
an external field) it is a consequence of Heisenberg uncer- 
tainty relations that the movement in this direction will 
be faster. The slow enough movement in the other direc- 
tions then will make up an adiabatic perturbation that is 
known not to change the state of the particle describing 
the fast motion. 

As a result, the wave function can be combined as a 
multiplication of an envelope depending only on the unre- 
stricted coordinates, and of the fast motion state depend- 
ing on the unrestricted coordinates as on parameters. 

This approach gives (see Appendix ^) for the eigen- 
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FIG. 2: Comparison of the first several spin wave modes 
frequencies in a spherical cavity obtained numerically (solid 
curves, JIt]]) and in the adiabatic approximation (dashed 
curves, Eq. (jss|)). The eigenfrequencies are plotted as func- 
tions of the ratio of the characteristic wavelength to the 
sphere radius. 



frequencies 



lu l + RVoj L + £Vw L 



(58) 

where u>l is the Larmor frequency in the center of the 
sphere, and n z ,n p = 0, 1,2, ...,oo are the longitudinal 
and radial quantum numbers respectively. Since £/R <C 
1 ~ a' n we see that the lower-lying levels belong to ao and 
therefore decay exponentially with diminishing z. Modes 
with n z = 1 will make a single oscillation before vanish- 
ing, modes with n z = 2 a double, etc. 

Eq. d58| ) reduces to (|50|) in the limit £/R — ► as it 
must. Indeed, in both Q58|) and ( |50| ) there figures the 
Larmor frequency at the top of a sample — z = L for 
(pOh and z = R for (58) and the second term in brackets 
in ( p8[) tends to zero when — > 0. 

If we introduce dimensionless frequencies according to 
jlTj u> a = u>l + RVuJLfa.) we conclude that the observ- 
able eigenfrequencies for our problem are described in the 
adiabatic approximation by 



/a = l 



i 

R 



(59) 



The comparison between this formula and the results ob- 
tained numerically in Ref. [|l7j is shown in Fig. ||. In 
Ref. [|l7j the combination £/R was designated as A. 
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FIG. 3: The dependence of the frequencies, half-widths and 
heights of the first four modes (numbered in the same order as 
on Fig. 2) on the strength of the demagnetizing field lom (solid 
curves) obtained by numerical simulations in Section V. Input 
parameters are = 0.1, C = 20. First order perturbational 
corrections (dashed lines) calculated in the text belong to the 
region u)m /2RVlol <C 1. Perturbational calculations were 
done only for the modes with radial quantum numbers n p = 0, 
among the four depicted only the first and the fourth are of 
the type. Heights of the modes are impossible to calculate 
perturbationally. 



modes. 

Quite analogously to the case of a finite cylinder, we 
obtain for the corrections to the modes (see Appendix 



X , .sphere 



U>M 




75 • (60) 



And correspondingly for the corrections to the imaginary 
parts of the modes 



e sphere 



(61) 



Numerical constants <&„ 2 are the same as in (pq). 

The parameter determining the relative value of the 
dipolar field effects in a sphere is 



(62) 




In the next section we will solve the eigenvalue problem 
in a sphere in the presence of a dipolar field of an arbi- 
trary strength. It is interesting to compare the results of 
numerical simulations with the first order perturbational 
corrections written above. 

To this end the positions 5R(w Q — u>l) and half- widths 
Qu> Q of the first four modes obtained from both numer- 
ics and analytics are plotted on Fig. || for several U)m- 
Heights are plotted only as obtained numerically. As one 
can see from Fig. ||, perturbational approximation is sat- 
isfactory for the values of ujm /2R\7lol up to ~ 0.1, or, 
for values of ( |6^ ) up to ~ 0.4. The small discrepancy 
between numerical and analytical results even in this re- 
gion is accounted for by the restrictions of the adiabatic 
approximation used to fulfill analytical calculations. 

The spectra themselves calculated numerically for sev- 
eral values oiuiM are shown on Fig. ||. Modes weights are 
redistributing between adjacent modes with ujm growing. 
Not only the weights but also half-widths and positions 
of the modes change. 



Note that the numerical scheme developed in Ref. |l7j] 
requires more and more computational effort for £ tend- 
ing to zero. The calculation time to get safe eigenfre- 
quencies values grows. That is why the numerical curves 
are not shown in the vicinity of zero. Since we ourselves 
use a similar computational technique we put off more 
detailed discussion until Section V. 

Contrariwise, the discrepancy between approximate 
and numerical curves at large £ is accounted for by inap- 
plicability of adiabatics out of the region £ -C R- 

So approximate and numerical methods complement 
each other. While numerics is the method of choice for 
relatively large £/R when adiabatics breaks down, it re- 
quires increasingly larger basis to obtain reliable results 
for small £/R. In this region it is easier to calculate 
eigenfrequencies in adiabatic approximation. 

Let us now look on dipolar field correction to the 



V. SPHERE. NUMERICAL CALCULATIONS 

The complexity of the Hamiltonian ( flc| ) does not allow 
finding exact eigenfunctions by analytical methods apart 
from using perturbation theory. Nevertheless it is always 
possible to solve equation ( |43| ) numerically. We chose a 
spherical container for numerical investigations. 

Before proceeding to the description of the simulation 
scheme utilized, a formalism is to be established for solv- 
ing an eigenvalue problem with a non-Hermitian Hamil- 
tonian. Of interest for us is the generalization of the ex- 
pansion ( p^ ) of the Green function into series over eigen- 
functions. The following subsection is dedicated to the 
topic. 

In the two remaining subsections we discuss respec- 
tively the technique and the results, other than already 
mentioned in the previous section, of the numerical sim- 
ulations. 
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FIG. 4: Spin wave spectra obtained numerically in Section V 
for several ujm /2RVujl < 1. Input parameters are £/R = 0.1, 
C = 20. Modes are numbered in the same order as on Fig. 2. 



It turns out then that, notwithstanding that neither 
the set | a) nor \a) is orthogonal, there is orthogonal- 
ity between the two sets. Indeed, following a conven- 
tional scheme of proving mutual orthogonality of eigen- 
functions, we note that 

(a\H\P)=Lu p {a\p). 

On the other hand, 



{a\H = [n + \a) 



(w*|a)) + = w a (a|. 



(65) 



Multiplying this by \(3) and subtracting the previous re- 
sult, we see that 



J* - up){a[ 



0. 



(66) 



Thus and \a) are orthogonal if a ^ ft. So it is said 
that the two sets |a) and \a) constitute a biorthogonal 
set of eigenfunctions. 

The expansion of an arbitrary function into a conver- 
gent series is then possible 



(67) 



Here the eigenfunctions are supposed to be normalized 
so that 



(a\ 



la/3- 



Performing such an expansion for \h + ) in 
obtains 



E 



\a)(a\ 



(68) 
]), one 

(69) 



Green function of a non-Hermitian Hamiltonian 



In this subsection we derive an analogue of (|45|) for 
a non-Hermitian Hamiltonian H ^ H + . For H ^ TL + 
the set of eigenfunctions \a) of 7i is not orthogonal. A 
second set of functions, viz., the set of eigenfunctions of 
the Hermitian conjugate operator Ti. + , is to be introduced 
& 

The eigenvalues of the operator 7i + are just the com- 
plex conjugates of the eigenvalues of 7i. Indeed, the 
eigenvalues are found from a secular equation. And for 
the Hermitian conjugate operators these algebraic equa- 
tions may be shown to be complex conjugate. 

Thus the set of eigenfunctions of Ti + may always be 
numbered with the same index a. To distinguish this set 
from | a) we will denote it by |a) 



H\a) 
H + \a) 



(63) 
(64) 



However, it should be borne in mind that \a) like \a) is an 
ordinary set of ket-vectors, which has the corresponding 
set of bra-vectors. 



This is the sought-for generalization of fl45j). For a Her- 
mitian Hamiltonian the two sets coincide |a) — \a}, and 
reduces to (^5|). 
Thus we see from ( fi"l"|) that the spectrum of a non- 
Hermitian Hamiltonian consists of a set of Lorentzians 
at oj = ^Ruj a with half-widths 3(jJ q each entering with a 
weight (h e+ \a){a\h e+ ). If the rf field can be regarded as 
uniform on the scales of the sample, the relative weights 
of the Lorentzian peaks are 



;r|a)d 3 r / (a|r')dV 



(70) 



In our particular case, the Hamiltonian (40) is symmet- 
ric T-L + = H* . This is trivial to see without the dipolar 
field term, but this term can also be shown to be (real) 
symmetric, since partial differential operator d\ and the 
Green operator Goo are both real symmetric and com- 
mutative d^Goo — Goodl (We remind that for an integral 
operator with a kernel G(r, r') Hermitian conjugate has 
the kernel G*(r',r)). Indeed, we write 



#&om = # / ^Xd 3 



m(r'R 2 



1 



(71) 
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Taking the integral two times by parts and taking into 
account that m(r')/|r — r'| — > at z' — > ±00, we get 



c^,m(r') 



dV = Qoodim.. 



(72) 



Therefore from the definition (|4j) we immediately 
conclude that in the case of a symmetric Hamiltonian 
H + = H* the two sets of eigenfunctions are related 
through (r|a) = (r|a)* = (a|r). Expression (70) for 
the relative weights of the modes simplifies then to 



'r\a)d r 



(73) 



Note that the (complex) value itself of the integral is 
squared, not its absolute value, as it would be should the 
Hamiltonian be Hermitian. We will use the expression 
( [73]) for the modes weights in numerical calculations of 
the spectra. 



B. Numerical approach 

One of the methods for solving a spectral Sturm- 
Liouville problem (conceptually, perhaps, the simplest) 
consists in its finite-dimensional approximation. For- 
mally, we then are left with standard algebraic spectral 
problem. For the solution of the latter one can implement 
one of ready safe well-established algorithms. However, 
one must be cautious with the dimension of approxima- 
tion. 

One of the possible discretization techniques is to find 
eigenfunctions in the representation of some complete or- 
thonormal set of functions when the Hamiltonian jiOj 
would become a matrix. Such a scheme was developed 
in application to spin waves in Ref. jl^], JT7] . 

A handy orthonormal set to choose is that of eigen- 
functions of the Laplace operator satisfying the boundary 
conditions for the geometry in question 



= 0, 



[d 2 + kl\ |m) 

rkd%\n)\ a = 0. 



(74) 
(75) 



Here /i stands for a complete set of indices needed to 
describe a state, are waveconstants. In the case of 
such a choice of the set the first term of the Hamiltonian 
( [Io|) becomes trivial and the boundary conditions are met 
automatically. 

The eigenfunctions |a), |a) of the operators H, H + 
then take the form 



*)=X»Ha>. 



(76) 



where the coefficients of expansion (fi\a) and (v\ca) are 
found numerically as right and left eigenvectors of the 
matrix form of the Hamiltonian ( ]40| ) corresponding to 
eigenfrequencies u> a 



^W^Ha) = w a (fi\a), 

V 

y]{a\v)Huti = w a {at\(J,), 



(77) 
(78) 



where = (fi\H\u). 

For a sphere /i denotes the set n, I, m of the radial, 
polar and azimuthal quantum numbers: n, I = 0, 1, 2, 
. . ., m = —I, + The corresponding basis 

is 



(r\n) = (r\nlm) = Cnij^kni^Y^fj) 



(79) 



where ji is the spherical Bessel function, Y™ is the spher- 
ical harmonic, and renormalization coefficients c n i are 
defined according to 



CniCn'i / ji(k n ir)ji{k n >ir)r 2 dr = 5„ 



(80) 



The waveconstants fc n ; depend on the boundary con- 
dition. That of (|75] ) requires that fc„/i? be the (n+ l)-th 
zero of the derivative d r ji (r) of the spherical Bessel func- 
tion. 

Since we are interested only in axisymmetric modes, 
which couple to a homogeneous rf field, we may simplify 
the formulas by working with a sub-basis \nl0). 

Infinite indexing n, I is to be truncated at some finite 
values for numerical computation. Maximum values of 
"max, ^max are restricted by computational tractability 
of resulting matrices. 

On the other hand, justification for such a truncation is 
that coefficients (nl\a) tend to zero for large n, I because 
of the oscillating character of ji. We expect (nl\a) close 
to zero when the characteristic scale £ of the function 
(rja) becomes greater then the period ~ R/n of the os- 
cillations of the basis function ji. Empirically, n max = 10 
is already quite good for customary £/R ~ 0.1. Note 
that the change in ji(k n ir) with increasing I is much less 
dramatic. So more I's are to be retained in the sub-basis. 
We used l max = 51. Further increase of n max , Z max proved 
to have no apparent effect on the spectra for £/R ~ 0.1. 

However, for £/R — > more and more sub-basis func- 
tions should be kept, which leads to rapid slowing down 
of the computations. In this limit adiabatic approxima- 
tion (see Sec. IV) gives safer results. 

It is convenient to seek for the eigenfrequencies ui a in 
the form ujl + RVu>Lfa- The Hamiltonian for the matrix 
equation on f a 



il\H f \n'l' 



'1\ 3 



+ (nl 



2RVul 



l + £) {k nl RY5 nn ,8 w 
n'l') (81) 
(nl J ^ — ri zz J TiT) 



comprises three parameters — £,/R, the ratio 
<x>M/2i?Vu>L of u>m to the total field gradient over 



the sample and the regime parameter C . Here we have 
used 1 — 2ri^ ere = |. 

The matrix elements of z/R and of n zz by integrating 
over the solid angle with appropriate spherical functions 
reduce to integrals over the radial coordinate r. 

The matrix elements of z/R were written in Jl7| , the 
integrals arising should be calculated numerically. The 
rather lengthy calculations of the matrix elements of fi zz 
were separated into Appendix [e| 

We cite here only the results. The matrix elements of 
z/R = rcos8/R, where 9 is the spherical polar angle are 
non-zero only if I' = I ± 1 



(ni| z/R \n',l' = l±l) 
f R 

= C x CniC n n> I ji(k n ir)ji>(k n >i'r)r dr/R 
Jo 



(82) 



Here A is the greater of I, V and 
\/\f±>? 



(83) 



The matrix elem ents of n zz are non-zero only if V = 
{1,1 ±2} (see ( p28| ) in Appendix g) 



(n l\\-n zz \n'l) = S n 



(«?) 



-t+V 



(84) 



(n l\\-n zz \n',l' = l±2) = c°c°_iC„;c„^i? 2 

k n iji+i(k n iR)ji(k n 'i'R) - k n >i>ji+i{k n >i' R)ji(k n iR) 
X k 2 - k 2 

K nl K n'V 

It may be verified that the matrix (nl\Hf\n'V) is indeed 
symmetric. 

The algebraic eigenvalue problem for ( |si| ) with (g^) 
was solved using standard subroutine from the Linear 
Algebra Package LAPACK. Eigenvectors were then nor- 
malized and the left and right eigenvectors used to find 
the modes weights. 

In calculating modes weights using (]77|), (|7^ ) we note 
that (r|000) = 1/W, where V = 4nR 3 /3 is the sphere 
volume. Thus 

J (r\nlm)d 3 r = W(000\nlm) = WS n0 S l0 S m0 . (85) 



The general formula ( |70| ) in this case reduces to 

y(000|a)(a|000). (86) 

C. Results of simulations 

Numerical results in the absence of the dipolar field 
were obtained in Ref. [0. When ujm = there remain 
only two parameters in the problem — the ratio of 
the characteristic wavelength to the sphere radius and the 
regime parameter C which determines half- widths of the 
modes. The dependence of the spin wave spectrum on 
in the limit of weak demagnetizing field was plotted 
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FIG. 5: Absorption signals of a polarized Fermi liquid in a 
spherical container for several regime parameters C in neglect 
of the dipolar field effects ujm = (grey filled curves) and for 
a dipolar field of an intermediate strength ujm /2-RVwz, = 1 
(solid curves). The ratio £/R was put 0.1. 



in Fig. H to compare adiabatic approximation with sim- 
ulations. A typical absorption signal in the absence of 
the dipolar field is depicted as grey filled curves in Fig. || 
for several C. Decrease in C results in modes broadening 
without changing their positions. 

When ujm is small enough the dipolar field constitutes 
a perturbation to conventional Silin spin waves. It is this 
regime which was studied perturbationally in Sec. IV. 2. 

In this subsection we shortly consider other results of 
numerical calculations, viz., regimes of intermediate and 
strong demagnetizing fields. Although these regimes were 
not realized so far in Fermi liquids (see Table I), one 
cannot leave out what is to be expected. 

As u>m increases the spin wave spectrum undergoes 
crossover from Silin type for small demagnetizing fields 
(ujm /2R\7uj l < 1) to magnetostatic type for large de- 
magnetizing fields (% /2RVuJl > !)• On a gross scale 
this transition is represented in Fig. [gJ for £/U = 0.1 
and C = 20 (grey filled curves) and C = 5 (solid curves) . 
Modes weights change so that adjacent Silin modes group 
into fewer magnetostatic modes. These separate at even 



larger ljm until a uniform (Kittel) mode singles out at 
extremely large u>m- 

This latter mode is the only one to remain because 
we chose a uniform radio-frequency field for the response 
of the system. And a non-uniform rf field is required 
to couple to non-uniform magnetostatic modes since the 
influence of external field gradient Vu>l is negligible for 
large u M . 

The behavior described is not altered by larger dissi- 
pation (smaller C) other than Silin modes grouping be- 
comes more pronounced (see solid curves on Fig. ^). 

The dependence on C of a spectrum for a demagnetiz- 
ing field of an intermediate strength (wm/2-RWl = 1) 
is plotted in Fig. || as solid curves. Pronounced adjacent 
modes for larger C merge into magnetostatic conglom- 
erates with no distinction for smaller C. No apparent 
relation, especially for smaller C, can be seen with the 
spectrum in the absence of the demagnetizing field. 



VI. DISCUSSIONS 
A. The method 

We studied the influence of the dipolar (or demagnetiz- 
ing) field on the spectrum of linear standing spin waves 
in a polarized Fermi liquid in a finite container. 

A somewhat resembling problem was studied in the 70s 
in ferrimagnets. deWames and Wolfram |25| considered 
not the Larmor-precession case ((]]) , but the true Landau- 
Lifshitz equations of inhomogeneous magnetization dy- 
namics that contain an additional exchange torque term 
depending on the Laplacian of the magnetization. This 
term is similar to that arising from the Leggett equations 
(see Eq. @). 

The situation in ferrimagnets and in polarized Fermi- 
liquids belong however to different limiting cases. In ferri- 
magnets the dipolar field effects — magnetostatic waves 
— are pronounced and exchange is a perturbation. In 
a paramagnetic Fermi liquid under conventional experi- 
mental conditions, on the contrary, it is the dipolar field 
whose effect is smaller. 

There is another obstacle of transplanting the ap- 
proach of deWames and Wolfram to our needs. In order 
to accommodate all the inhomogeneities — both from the 
Landau-Lifshitz equations and from the Maxwell equa- 
tions (§) — the resulting differential equation for magne- 
tostatic potential is of the sixth order in spatial deriva- 
tives. Though it was manageable in a quasi-lD situa- 
tion (ferrimagnetic substance yttrium iron garnet YIG is 
grown as thin films), it becomes virtually intractable in 
any, even the simplest, 3D geometry. 

An approach to include dipolar field into magnetiza- 
tion dynamics analogous to ours was used by Deville et 
al. |g9fl in application to solid bcc 3 He. In this mate- 
rial the dipolar field gives rise to multiple spin echos at 
times nr following two isolated rf pulses at times and r. 
Ref. p9| explained this phenomenon quantitatively using 
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FIG. 6: From bottom to top — evolution of the spin waves 
spectrum of a polarized Fermi liquid at a crossover from weak 
to strong demagnetizing field. Several adjacent Silin modes 
form a group which evolves further as a magnetostatic whole. 
The distance between magnetostatic modes grows with wm 
until we are left with the sole mode in which magnetization 
in the whole sample oscillates uniformly. Input parameters 
were £/R = 0.1, C = 20 (grey filled curves) and C — 5 (solid 
curves) . 



a simplified local approximation for the dipolar field valid 
for a slab infinite in the directions perpendicular to the 
external field (although this restriction was not properly 
emphasized in the paper). Later Fomin and Vermeulen 
ppf utilized the form of the dipolar field term of Ref. 
to study dipolar corrections to a two-domain coherently 
precessing structure. 

As for linear spin waves, in this local approximation 
neither the relative positions of the modes nor their half- 
widths and heights change, and the spectrum only shifts 
as a whole. The demagnetizing field is by definition 
shape-dependent and it is erroneous to use the local ap- 
proximation a priori. 
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TABLE I: Comparison of the conditions of different spin-wave experiments Jio| |, Jig], JT^ |, [ |l8| in pure 3 He and in solutions 
of 3 He in 4 He. Prom left to right, i) The concentration x of 3 He atoms in mixtures (for pure 3 He not given), ii) The static 
magnetization M in the units of the equilibrium magnetization Mo (^) in the external field; different from 1 only in jlo| |, | ]18| |; 
see the text for more explanation, iii) The Larmor frequency uil/^tt, where ujl = r yH e . iv) The gradient X7H e = Vi^l/i of 
external field; v) The mean distance between modes £\7u>l in the units of frequency, £ being the characteristic spatial distance 
( [l9| ) — the wavelength of Airy type spin wave, vi) lum = AtyjM — a characteristic of the static demagnetizing field, vii) The 
radius R, for a box this is half-size, for factually used cylinders it is also roughly the half-height, viii) - x) The three parameters 
entering the Hamiltonian dil] ) used for simulations of the eigenvalue problem in a spherical cell, (£/-R) 3 and ojm /2RSJujl give 
the relative importance correspondingly of the exchange molecular field and of the demagnetizing field to the gradient of the 
external field; in a box these parameters lack strict meaning and are presented here in parentheses only for an estimate; C 
is the regime parameter (jij) at the temperature shown for each experiment in the left column, xi) The dipolar parameter 
( jB7j ) (cylinder and box) or ( [jj8| ) (sphere) (more exactly, their more accurate equivalents (^) and ( |6C| ) respectively); if it is 
much smaller than unity it signifies that perturbation theory can be applied at the conditions of the experiment. The input 
parameters for the calculations come from various sources fig]. 
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B. Main results 

The strength of the demagnetizing field is propor- 
tional to the polarization of the liquid. So in conven- 
tional weakly-polarized liquids the dipolar field can be 
neglected, while for strong enough polarizations its ef- 
fects overwhelm usual Fermi liquid exchange field spin 
dynamics. 



Because of its long-range character, the dipolar inter- 
action introduces an additional non-local term into the 
equations of the spin dynamics. This term is an integro- 
differential operator on the magnetization, wherein the 
integration is taken over the volume of the liquid. As a 
consequence the corrections to the spin dynamics depend 
strongly on the shape of the container. 



In particular, the influence of the demagnetizing field 
on the spectrum of standing spin waves in an (infi- 



nite) slab reduces for an arbitrary strength of the dipo- 
lar field simply to uniform shift of all the modes by 
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—u>m — —AirjM, where M is the static magnetization 
of the liquid in an external field, 7 — gyromagnetic ra- 
tio. 

In a finite volume of a liquid strengthening of the de- 
magnetizing field results in the crossover from the Silin 
type spin wave regime to the regime of magnetostatic 
waves. In this latter the mechanism of the forming of 
the standing spin waves has nothing to do with the ex- 
change physics of Fermi liquids. Magnetostatic spectra 
are in no way specific to the Fermi liquid and so contain 
no information on its parameters. 

For intermediate polarizations the demagnetizing field 
in a finite volume of a liquid leads to the corrections to 
Silin type spin wave spectra, both shifting the spectra 
uniformly and also changing the distances between the 
modes and the modes widths. These changes of the sec- 
ond type which distort the spectra are undesirable from 
the point of view of deriving parameters of the liquid 
from the spectra. 

We have seen that the parameter determining the in- 
fluence of the demagnetizing field on the spectra of spin 
waves in a finite cylinder in the first order of perturbation 
theory is 



2w 



M 



7t£Vw£, R 



log- 



(87) 



where ujm = 4:irjM characterizes the magnetization den- 
sity, Vujl is the gradient of the Larmor frequency, R is 
the radius of the cylinder base, L its height and £ is the 
wavelength ( fl9| ) of an Airy-type standing spin wave. The 
quantity ^Vwl gives the average distance between modes 
in the units of frequency. <f> and \& are numbers of the 
order of unity. 

For a sphere, an analogous parameter was 



The values of these parameters (along with the values 
of some others having appeared in the text) for several 
recent spin-wave experiments are arranged into Table Q. 

The ratio M/Mq of the absolute value M of the static 
magnetization to the equilibrium magnetization Mq ac- 
counts for the possible higher then equilibrium polariza- 
tion of the liquid. Polarizing a liquid M/M times its 
equilibrium value increases proportionally the strength 
of the demagnetizing field as well as the parameters ( |87| ) , 
(||). For the experiments @, where M/Mq could 
be changed, we present in Table I the data for several 
integer values of M/Mq which are close to the real ex- 
perimental values. 

For conditions when dipolar parameter (|^) or ( |8g| ) is 
smaller than unity one can use perturbation theory to 
get corrections to the modes frequencies. 

As the comparison between perturbation theory based 
analytical calculations and numerical simulations for a 
sphere shows (see Fig. ||) perturbation theory has ac- 
ceptable accuracy up to the values of to m / '2i?V 'to % ~ 0.1 



or (p7j), (|88T) ~ 0.4. However, at greater values of the pa- 
rameters^^) 1 (H) the discrepancy between perturbation 
theory and numerics grows dramatically. 

When the spin wave spectra are used for measuring 
the transverse relaxation time t, a proper treatment of 
experimental data taking into account dipolar field cor- 
rections is necessary. We estimated that the dipolar re- 
strictions on the correct determination of r from the con- 
ventional interpretation of the spectra are of the order of 
the parameter ( |8^ ) or (]8q). In particular, for the experi- 
ment |l(J about 4.2%. The latter means that the effect of 
a dipolar field can not significantly change the main con- 
clusion of this article that the polarization induced zero 
temperature spin wave damping does not exist, which is 
in disagreement with previous spin echo experiments [||, 
§, |, §• 

A major inference for planning future experiments is 
the proposal to use ellipsoidal, in particular, spherical 
containers, not only because the estimation of the shape- 
dependent dipolar field effects is simpler, but also because 
there are two roughly equal contributions to the spin 
wave spectrum distortion, — one from the inhomogeneity 
of the static demagnetizing field and the other from the 
demagnetizing field produced by the rotating part of the 
magnetization. The first contribution exists only in non- 
ellipsoidal samples, in which a homogeneous static mag- 
netization produces inhomogeneous demagnetizing field. 
For this reason implementation of such shapes for the ex- 
periments on the elicitation of the liquids characteristics 
from the spectra is disadvantageous. 

The dependencies of (|87|), ( p8| ) on experimentally con- 
trollable parameters are as follows (p7|)oc lu l log ujl, oc 



log(D /«ri), oc Vw^logVwz,, oc R- 1 , and || 



OC LO 



5/4 



oc (Do/kti)- 1 / 4 , oc Vlj l ' , oc i?" 1 / 4 . So a bigger R 
and a bigger gradient diminish the contribution of the 
demagnetizing field. 

But for typical experimental conditions as one can 
see from Table ffl, strongly-polarized 3 He- 4 He solutions 
never are too far beyond the regime of Silin spin waves 
perturbed by demagnetizing field, whereas pure 3 He at 
strong polarizations is in the "deep intermediate" regime, 
for which the results of Sec. V apply. In view of this, 
even at most favorable cell size and field gradient, pure 
highly-polarized 3 He seems to be unsuitable for a study 
of Silin waves spectra. 
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APPENDIX A: SPIN WAVES IN AN INFINITE 
MEDIUM IN A UNIFORM MAGNETIC FIELD 

The non-local dipolar term n[m] is known to become 
local in the important case of a medium infinite in two 
directions, with the proviso that m depends on only the 
remaining third coordinate. The direction in which m 
varies we denote as s, and then the medium should be 
infinite in the two directions perpendicular to s. 

In such conditions 

didj J y^lrf*' = SiSjd. J d s g ld (s - s')f(s')ds', 

V 

(Al) 

where by definition 



d a 



d 2 r' 



Vr'x + a 2 



dy 



3/2 



-27rsgna 



Here r_L denotes the coordinate vector in the plane per- 
pendicular to s. 

So n.ij[f] = SiSjf and 



n[m] = s(sm). 



(A2) 



The case of an infinite medium in a uniform magnetic 
field (i.e. Vwl = 0) is the simplest. Looking for a solu- 
tion of ( [si] ) in the form of a running wave 



m(r) = m e 



ikr 



(A3) 



we have s = k. Hence, Eq. ( p4| ) from which the com- 
ponents of the constant mo are to be found, becomes a 
linear algebraic equation 



w L + u M k'l -iut + uj M k x k y 
iuj + Lo M k x k y Z) L + u M k 2 



m 



47T ' 



where ujl = — Dk 2 and uj = to + i(D/C)k 2 . 

mo as a function of frequency has resonance at the uj 
which render the determinant of the matrix zero. This 
gives a Holstein-Primakoff type spectrum with an addi- 
tional attenuation term due to dissipation C _1 ^ 

id = \J (uj l - Dk 2 ) (uj l - Dk 2 + lo m sin 2 9) - i^k 2 , 



where 9 is the angle between z and k. 



(A4) 



APPENDIX B: SPIN WAVE SPECTRUM 

A specimen placed into the field of an NMR coil 
changes its impedance in two ways. First, the induc- 
tance L alters because so does the average energy of the 
field || 

^r = ^f H(t,r)B(t,r)d 3 r (Bl) 



due to dispersion. / is the current through the coil. A 
line over an expression here designates time average over 
oscillation period. 

Secondly, the resistance R appears owing to the dissi- 
pation of the energy of the field in the specimen 

RI 2 = -3- / H(i,r)5 t B(t, r)d 3 r. (B2) 

47T J 

Making use of the general solution (p|)-(pT|) and of the 
expansion (|24|) we write 



H(t) = H e z + h e (t)+H dip , (B3) 
B(t) = (i7 e + 47rA/)z + h e (i)+47rm(i)+H dip (B4) 

In taking space integrals of mutual scalar products of dif- 
ferent terms of H(£) and B(i) we note that those contain- 
ing Hdip transform into integrals over a remote surface. 
Since beyond the specimen magnetization is zero, such 
integrals vanish. 

In bilinear expressions we should write the monochro- 
matic rf field h e (t) = hf n coswi as h e (t) = i(h e + h e *), 
where h e = h^e _4Wt . Similarly, the rotating part of 
the magnetization should be written in the form m(i) = 
i(m + m*), where m = m m c-'( ui+ ^(l,i). 

Then the impedance Z = R — iu>L/c 2 of an NMR coil 
can be written in the form 

Z = Z a +r 2 J h e (t, r) (d t -iuj)m(t,r)d 3 r 



Z - iujr 2 / h e *md 3 r 



(B5) 



where Z is the impedance without sample. So the 
change in the impedance of the coil due to sample is 
proportional to 



2 J h e *md 3 r « (h 



e+l m+) 



(B6) 



The real part of this quantity gives the dispersion spec- 
trum, while imaginary — absorption. 

Introducing Green operator fl42| ) and normalizing, 
we arrive at expression (|4l|). 



APPENDIX C: DIPOLAR CORRECTIONS TO 
MODES IN A FINITE CYLINDER 

In this Appendix we derive expression ( |56| ) for mode 
shifts due to dipolar field in a finite cylindrical cell. 

The normalization coefficients of transversely homoge- 
neous rip = 77i~0 modes ([55]) are 



irR 2 ^ 



L/S- 



Ai 2 (a 



dx. 



(CI) 



The upper limit may be put equal to infinity and then the 
last dimensionless integral is a number depending only on 



n z . 
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In calculating the integrals in (|5J) we use the expansion 
of the Green function in cylindrical coordinates j22[ p. 
140] 

= y- e imW-<p') / e -*k-«'| j m (k p )J m (kp')dk. 

(C2) 

Integrals over <£> and 99' give (27r) 2 <5 m o- Then the inte- 
grals over p and p' with J give R 2 Jf(kR)/k 2 so that 



/ * (z)n zz Mz)]d 3 r 
Jv 



00 L 
2C 



- 7tR z 



J({kR) 
k 2 



$(z) / $(z')afe- fc|z - z ' l dz'^fc. 



In the first integral in @ = ipl^(z), $(z) = 1, 

while in the second — ^(z) = &{z) = ^n-{z)- 

As a result of differentiating expansion dC2|) we have 



(C4) 



d 2 e -k\z-z'\ = _ 2k s( z - z ') + k 2 e- klz - z ' 1 . 



The integral of the ^-functional part is the simpler, using 
/ x~ 1 Jl(x)dx = i we obtain ttR 2 <S>{z)$(z)dz. For 
both integrals in ^5J) this gives unity. 

So the (5-functional part in the dipolar operator for 
transversely homogeneous spatial distributions gives the 
local slab value (p2[). 



The second part in (C4) is shown below to be non-zero 
only for finite samples. Indeed, it yields 



-irR / / $(z')F 

Jo Jo 



R 



dz'dz, 



TABLE II: Numerical constants in the expression ( p6[ ) for 
mode shifts due to dipolar field as functions of the mode num- 
ber. 
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And for the second integral in (|54| ) 



l + ^ilog^i 
7T R 5 8R 

Here the n z -dependent constants 
Ai(x)dx) 2 



(C8) 



(C9) 



log B„ a 



log <d, lz 



2 Ai 2 (x)dx ' 

J^°, xAi 2 (x)(ia; 

n? Ai 2 (x)dx ' 

I- a ' n x\og(x)Ai 2 (x)dx 
f^°, xAi 2 (x)dx 



(Cll) 
(C12) 



where 



are of the order unity, as is seen from Table ||, where 
they are calculated numerically for the first six modes. 

Plugging (C7) and (C8) into we get for the dipolar 
corrections to_the modes frequencies in a finite cylinder 
expression 



F{p) = / J?(x)e- px dx 
Jo 



(C5) 



is the Laplace transform of J 2 . Although its value can 
be found in tables (see, e.g. p5| , formula 6.612]) for our 
purposes it is sufficient to know its value for small p, 
where it diverges logarithmically 



F 



R 



1 e 2 \z-z'\ 
log ■ 

7T 



8R 



(C6) 



Writing for small z 



L 1 e2 \ z — z '\ , 1 r , e L , z m ?\ 
log — — -dz 1 = L log — + z log — + 0(z 2 ) 



n 



8R 



8R 



eL 



and substituting ^j 2 - — > x we obtain for the first integral 



in (54) 



1. Dipolar error to transverse relaxation time 

We are now in a position to estimate the error in the 
determination of the transverse relaxation time because 
of the dipolar field. 

As an example, we consider the experiment |fio|| , 
wherein r was obtained from the regime parameter C, 
which, in its turn, was determined from the slope of the 
linear fit to the dependence of modes half- widths ^su a on 
their position 3?(cj q — u>l — LS7ujl)- In the collisionless 
regime in a finite cylinder, as is seen from j50| ) and (f3ll), 
the two quantities are related through 



1 

3C' 



(C13) 



We now find the dipolar correction to this value. 

For finite C _1 the scale £ should be replaced with 
the complex £(1 + i/C) 1 / 3 . This means that the dipolar 
field changes also imaginary parts of the eigenfrequencies 
and hence half-widths of the modes. In the collisionless 
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regime (C> 1) the imaginary part of the correction (|5 
equals 



8i? 



(C14) 

So the dipolar field changes the widths of the modes 
proportionally to ojm- Since for all parameters being pos- 
itive the initial imaginary parts ^aj, ^Vwi of the modes 
( po| ) are negative, the modes narrow down in the first ap- 
proximation. 

As a matter of fact, experimentally measured are not 
the absolute positions of the modes frequencies but rather 
their positions relative to each other. So modes positions 
in Ref. [0 were determined relatively to the position u> 
of the zeroth mode: not ui a but ui a — luq. 

To diminish the error due to the data scattering in the 
value of the slope derived from the fitting, it is desirable 
to fix the Larmor frequency at the wall wl +LVu>l which 
is experimentally badly determinable. So the position ujq 
was put equal to its value = LOL + L\7ujL+a' ^VujL in 
the absence of the dipolar field. This mean s tha t instead 
of uj a — lol — LVujl in the denominator of (C13) actually 
used in Ref. Jl(J was $l((u> a — ujq) + a^Vwt) — 

wo) -wl- LVuj l + 

Calculating the ratio of $su) a to this q uanti ty taking 
into account the dipolar corrections (56), (C14) we get 



3CV a' n £Vu L nR 



$ log 



8R 



*o log ■ 



e £ 



(C15) 



So the error introduced to the determination of C is of 



the order of (87). 



APPENDIX D: SOLUTIONS FOR A SPHERICAL 
CELL IN ADIABATIC APPROXIMATION 

In this Appendix we obtain spin wave solutions in a 
sphere neglecting the dipolar field. 

For a solution of the three-dimensional eigenvalue 
problem in a sphere of a radius R 

[Dd 2 + zVuj L ]u(z,r ± ) = 5uiu(z,r ± ), (Dl) 



d r u 



r=R 



0. 



(D2) 



where we denoted Suj = lo — u>l for brevity, the adiabatic 
approximation consists in the substitution 

u(z,r ± )=v(z;p)w(p)e im '<>, (D3) 

where v nz (z;p) are the eigenfunctions of the equation 

[Dd 2 z + zWuj l ] v„ z (z; p) = oj nz (p)v n , (z;p). (D4) 

The coordinate notations are represented in Fig. [?]. 
The appropriate boundary conditions for (D4) will be 
discussed below. 




FIG. 7: A scheme of the coordinate system used to solve the 
eigenvalue problem in the adiabatic approximation. 



Then Eq. ( |Dl| ) becomes 

v nz (z; p) [D&1 + UJ n% (p) - Suj] w(p)e m ^ 
(wd 2 ^v„ z +2d p wdpV nz ] 



0. (D5) 



Adiabatic approximation utilizes the fact that the sec- 
ond term could be for certain conditions neglected with 
respect to Dd 2 _w(p). Then we would obtain an effec- 
tively decoupled eigenvalue problem describing the slow 
transverse motion 



[DDI + U) n M} w{p)e mip = 5uw{p)e miv . 



(D6) 



For an estimate of the conditions of the possibility of 
that neglect we write 



d p Vn z {z;p) ~ d p zd z v Uz ~ (p/R)d z v nz . 



(D7) 



and from (D4) we see that d z v Uz ~ Vn z /£, where £ is the 
Airy spin wave characteristic wavelength (^9|). On the 
other hand d p w ~ w/p. So, the condition of applicability 
of the adiabatic approximation to this problem is that the 
ratio of the second term in ( |P5| ) to Dd 2 _w(p) should be 
smaller than unity 

p 2 /R£ < 1. (D8) 

The boundary condition ( |P^ ) on u reads 

[cos 9d z + sin dd p ]v nz (z- p)w{p)e lmifi = 0. (D9) 

For small we have cos# w 1 — p 2 /2R 2 and sin# ps p/R. 
Estimating d p v nz as above, we see that sin9d p v nz is of 
the order of (p/R) 2 d z v nz , and in first order in p/R the 
boundary conditions are 



wd z v n ^ + ^rv n ^dpw) 2 



0. 



(D10) 



The first term is proportional to l/£, while the second 
to 1/R and in the case £ <C R the boundary conditions 
reduce to a much simpler form 



<zxR- 







(Dll) 
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Finally we formulate once again all our assumptions, 
i.e. p, £, p 2 /£ <C R. These three reduce to £ <C R if 
p ~ where l/^/a is the characteris tic tr ansverse 

spatial scale of the wave function (see Eq. ( D16 ) below) . 

Now we proceed to the solution itself. First, Eq. (P4 ) 
should be solved with the boundary conditions (i) ( pll ) 
at z = R — p 2 /2R and (ii) — > as z — > — oo. 

The solution is then 



v nz {z;p) = Ai 



(p) 



where 



w„ = (p) 



R 



EL 

2R 



(D12) 



(D13) 



and a' n < is the n-th zero of the derivative of the Airy 
function Ai'. 

Then Eq. ( |P6| ) describes a two-dimensional harmonic 
oscillator 



Ddl 



Vu L p 2 



w(p)e 



imLp 



(D14) 



i? 2 

Solution in polar coordinates is 

w npm (p) = P lml e- a "^ 2 L^(ap 2 ), (D15) 

where L%(z) are the Laguerre polynomials and we de- 
noted 



9 Vwl 

a" = 



2DR 2£ 3 i?' 
The corresponding eigenfrequencies are 



(D16) 



lu l + RVlu l + £Vcj l 



R (U + 1) 



where wl is the Larmor frequency in the center of the 
sphere, n = 2n p + \m\. 

Again we remark that only the modes with zero az- 
imuthal quantum number m couple to the rf field. Eigen- 
functions with m — occur for even n = 2n p , as is writ- 
ten in 



1. Effects of dipolar field in the first order of 
perturbation theory 

We now calculate the demagnetizing field corrections 
( pi) ) to the modes in a sphere. 

As we saw in the previous subsection, the dipolar-free 
solution of the equations of motion satisfying boundary 
conditions (44) is 



^»,»„o = c nznp oAi( R -f- z + a'^)e- a P 2 / 2 L np (ap 2 ), 

(D17) 



where n z ,n p = 0, 1, 2, oo, and m was put equal to zero 
because only the modes with m = couple to the homo- 
geneous rf field. Furthermore, for simplicity we consider 
modes with n„ = for which L n = 1. 



In order to calculate the normalization coefficient c 



a , on 



Wi 2 00 



2tt / e~ 



1 L r> n 



Ai 



.2(R- 



-)dz pdp 



we make the substitution x = (R — p 2 /2R — z)/£ 



9 2-R 

oo z 5 



u n = 00 



2tt£ J e~ ap2 J Ai 2 (a: + a' n J dx pdp 



— I Ai 2 (x + aL ) dx, 



(D18) 



Here the upper limit of integration ov er x m ay be esti- 
mated as - (R/£ - l/aR£) = (R/£ - y/2j/R) > 1, after 
which the integral over x decouples from that over p in 
the approximation <C 1. 

In the case of a sphere n zz [1] = nS z ' = | and as was 
already mentioned there is no contribution to the spec- 
trum distortion from the static dipolar field. One needs 
only to calculate the second integral in 

Again, the S— functional term in (C4) using p|| 



kJJkp)JJkp')dk = -S(p - p') 
P 



can be easily seen to give a constant 1 as it should. In the 
limit £/R — > a sphere transforms into a half-space, the 
modes being localized near the boundary. A half-space 
is a particular case of a slab, with the height L — > oo. 
Hence the solutions in a sphere in the limit £/R — ► 
transform into solutions for a thick slab depending on 
only z. The demagnetizing tensor for such solutions re- 
duces to a constant (f>2). 



The second term in (|C4j) gives 



V 

oo oo oo 

x f e - ap ' 2/2 J (kp') [m(x + a' n ) [ki{x' + a' n ) 



x e 



2 _ 12 
X ~ X +P 2R( 



dx' dxp'dp' pdpdk. (D19) 



The two terms under the module sign appeared from \z- 



z'\. We can estimate x 



1 and (p 



- P ,2 )/m ~ 

x'\. Then the 



l/2R£,a = y/£/2R, therefore \z-z'\ «J% 
integrals over p and p' can be taken |2g, formula 6.631] 



/ e- ap2 / 2 J Q (kp)pdp=-t 
Jo a 



-k 2 /2a 



(D20) 
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The integral over k 



where r< (r> ) is the smaller (larger) of r and r' . Then 



-k 2 /a e -kp k 2 dk 



(D21) Goo\nlm) 



4tt 



where p = £ | x — a/ 1 , is a funct ion of p^/a ~ y£jR. There- 
fore, we can substitute (D21) with the zeroth term of its 
expansion in series with respect to P\/a, which is j^/ttcl 3 . 
Then 

J i ; (z)n zz [t/} nz (z)} d 3 r = l- yW; 



= 1 - 



2~R' 



(D22) 



where are the same numbers as in the case of a cylin- 
der. 

In the end we obtain ( |6C| ) for the corrections to the 
modes. 



APPENDIX E: MATRIX ELEMENTS OF THE 
DEMAGNETIZING OPERATOR 

In this Appendix we are going into the detail of cal- 
culation of the matrix elements of the dipolar integro- 
differential operator fi zz . We show that only those el- 
ements (nlm\n zz \n'l'm') are non-zero that are between 
the states with to = to' and I' = I, I ± 2. 

Before proceeding we remark on notations. We will 
write the integral operator in (nj) as 



SooM(r) 



(El) 



justification being that Goo is the Green operator for 
Laplace equation with the boundary condition of van- 
ishing at infinity. 

So our plan for this section is, first, to calculate 
Goo | rilm). Then, secondly, we calculate the result of act- 
ing of d 2 on an arbitrary function /(r) expanded in spher- 
ical harmonics 



/(r) = £/ ;m M>T(f) 



(E2) 



im 



By substituting Qoo\nlm) for / we eventually find the 
matrix elements themselves. 



1. Calculation of Goo \nlm) 

We find Goo\nlm) directly by integration, using a well- 
known formula from the theory of spherical harmonics 

I 4-7T 7*^ 



21 + 1 



(E3) 



rl+ T+l ^m)3i{kmr'y 2 dr>. 



The integrals at the right-hand side can be taken easily 
as follows. We notice that IM 



l ji+i{x) = -d x [x l ji(x)] , 



(E4) 



= 8 X [x l+1 j t (x)] . (E5) 
Taking integrals from the two sides we obtain 



(|) ji-i(y)dy = ji(x), (E6) 
(~\ Ji+i(y)dy = jt(x) - Ma). (E7) 



On plugging the above into (E3) and using another 
property of the spherical Bessel functions 

21 + 1 

(x) + ji-i (x) = ji(x) (E8) 



we get a sum of two terms 

Goo\nlm) = n n m9 \nlm) 



(E9) 



(kmR) 2 

- 2T+-lk^R Yl (r) (i?) 

This expression is inapplicable when k n iR = 0, which 
takes place for n = I = 0. In this particular case, inte- 
grating explicitly, we get 



£oo|000) = 471-000*2 T(r)jo(0) 



R 2 



2. Calculation of 8%f(r) 

Making use of the overt expressions of the basis vectors 
in the spherical coordinates 

r = x sin cost/? + y sin sin ip + zcos#, 

= xcos 6 cos ip + y cos 9 sin ip — z sin 6*, (E10) 

tp = — x sin ip + y cos ip, 



and of the nabla operator 



d = W r + 6-dg + tp—^—dv, (Ell) 

r r sm 



it can be verified that 



d z = cos 9d r + — sin 9 (e^L - e^l+j ■ (E12) 
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Here 1 = — ir x d is the angular momentum operator 
h = -id v , l± = e ±lv (±d e + i cot 6d v ), (E13) 

which has the well known effect on the spherical harmon- 
ics 

l z Y™ = m >7\ (E14) 
l + Yt m = y /(l-m)(l + m+l)Y l m+1 , (E15) 

(E16) 



l-Y^ = ^(l + m)(l-m+l)Y t 



m—1 



The product of each spherical harmonic with sin 9e tip oc 
Y^, cos (9 oc Y?, or sinfle"^ ex Ff 1 , is a sum p?| 



sin(9e- I¥, K 



mym-l imym-1 



(E17) 

smOe^Yl" = -a^Yfi? 1 + bJ^Yffi 1 , (E18) 
cos^™ = cT +l Y^ l+ cTY{2 x (E19) 

of the harmonics with the adjacent I and m multiplied 
each by a coefficient (which, in fact, are particular cases 
of the Clebsh-Gordan coefficients) 



l-m+l l-m+2 
21 + 1 21+3 ' 

l+rn l+rn—1 
21+1 2Z-1 ' 



(E20) 



m _ / l+rn l — m 
C l ' V 2i +! 2i— 1 ' 



Hence from (|E12|) 

^/ = E (- c ?+i^+iA + + cr^iA") /«m(r). (E21) 



Here the coefficients turn out to be the same as in (El£), 
and we introduced two differentiating operators 



-dr + l/r, 



Ci = d r + (l + l)/r. 



(E22) 



We may rewrite (E21) by shifting the summation in- 
dices as 

d z f(v) = Y / (d z f)i m (r)Yr(r), 



where 

(dzf)lm = —C%'£l-\fl-l,m + c T+l^l + lfl + l,m- 

Then repeating the procedure, we get 

af/(r)=^(a 2 2 /) im (r)^(r), 



(E23) 



with 



„m m r+ r+ f 
I H-l'~l-l L 'l-2J l -' 2 S n 

(cTf + (cr+lf] {d^Jlr, 

to TO fi— r— f 
c l+l c l+2 J ~l+l i -l+2Jl+2, , m- 



I in 



(E24) 



TABLE III: Operators Cf , L~[ 



f(r) 


r l 


l/r l+1 


3i( r ) 


ji (const r) 







(21 + l)//r 


ji+iir) 


const ji+i(constr) 




(2H 


-l)f/r 


Ji~i( r ) 


const _ i (const r) 



Here we used that C'^_ 1 £ t = C l+1 C~l = — (<9 2 ) ; , where 
(<9 2 ) ; is the /-th component of the Laplace operator 

(n=^ 2 +4-^- (E25) 

Changing the summation indices in each term again 
we arrive at 



®*f — Sim ( c ™l c "+2^+2(f)^rflA^ 



(E2C) 



+ 



+ C lC r^ 2 (?)A- lA"J flm 

3. Matrix elements 

Substituting Goo\nlm) from ( [Eg) into (|E26|) and using 
le i 
find 



the properties of £ ; summarized in Table [iTj, we 



1 

4tt 



tftGaolnlm) 



= ^ +l cT +2 c nl3l+2 {k nl r)Y^ 2 {v) 
- '(cT) 3 + (^i) 2 ]Nm) (E27) 

+ cttcrc nl Y^ 2 (r)(ji. 2 (k n ir) 

Ji-i{KiR) (r\i-A 



(21-1)'- 



k n iR 



(5) 



It is not hard to verify by straightforward differentiating 
that this expression holds also for n = I = m = 0. 

The matrix elements (n' V m 1 \d^Q ^nlm) are to be ob- 
tained from this expression by integrating with the com- 
plex conjugate of ( 79|) over the spherical volume of our 
sample. In doing so we see that only the elements with 
to' = m, V = I, l±2 are non-zero as we already mentioned 
in the text 

(nlm\ £ fiJ&oKZ' = l-2,m) 

= (cT-icT) l R cc'MkrMk'rydr; 



(nlm\ i dlG^n'lm) = [(cf) 2 + (tfttf 
(nlm\ ± <9 2 £ 00 |n',Z' = / + 2,m) (E28) 

= (c^cpV (j* ji{kr)ji(k'r)r % dr 
, Jl+1 (kR)j l+1 (k'R)\ 



-(2l + 3)R 6 



kR 



k'R 
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where we introduced shorthand notations k — k n i, k' = 

kn'l'j c — Cnli £ — C-n'l' • 

In the last term we used 

/ ji{kr)r l+2 dr/R l = R 3 j l+1 {kR)/kR (E29) 
Jo 

with k = k n[ . There is one exception, though, when 
I = n = — then fc„z = fcpo = and the division in 



the right-hand side of (E29) is undefined. Integrating 
straightforwardly, we get instead R 3 jo(0)/3. 

At the e nd, w e take integrals in the off-diagonal el- 
ements of ( E28 ) using the formulas below |p5| , formula 
5.54] 



ji{ ax)ji(bx)x 2 dx (E30) 

2 bji^ 1 (bx)ji(ax) ~ aji_ 1 (ax)ji(bx) 
~ X a 2 ~b 2 

2 -bji+i( bx )ji (ax) + aji+i(ax)ji (bx) 
~ X a 2 ~V 

2 bji +1 (bx)ji +2 (ax) - aji + i(ax)ji +2 {bx) 
~ X a 2 -b 2 

+(2l + 3)x s Max)jl+ l ibx) . 

ax bx 

Here a ^ b and we consecutively applied recurrent rela- 



tion (E8) first to ji-i and then to j\. 

Plugging the integration limits we see that above- 
diagonal (V = 1 + 2) elements equal corresponding below- 
diagonal (I' = I — 2) with the appropriate change I' <-> I. 
(For n = I = this is to be verified manually.) 

Finally, we get formula in the text for the dipolar 
part of the Hamiltonian. 



The case a = b in ( E3C ) is used in the calculation of 
the normalization coefficients c n i (|80| ) 



/ 



jf(ax)x 2 dx = — (jf(ax) - ji- 1 (ax)ji+i(ax)) . 



We get 



Cnl = 



[jf{k n iR) - ji-i(k n iR)ji+i(k„iR)] 



-1/2 



R 3 
2 

^-[j?(k nl R)+j 2 +1 (k nl R) 
21 + 1 \~ 1 / 2 

-j ~-jl{knlR)jl+l{knlR)\ 

k n iR 



when n + I > 1. For n = I = we have koo = and 
ji-i grows unlimitedly when its argument tends to zero. 



Then we integrate directly coo = {R Jo/3 



-1/2 
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